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We discuss a toy model for adiabatic quantum computation which displays some phenomenological 
properties expected in more realistic implementations. This model has two free parameters: the 
adiabatic evolution parameter s and the a parameter which emulates many-variables constrains 
in the classical computational problem. The proposed model presents, in the s — a plane, a line 
of first order quantum phase transition that ends at a second order point. The relation between 
computation complexity and the occurrence of quantum phase transitions is discussed. We analyze 
the behavior of the ground and first excited states near the quantum phase transition, the gap and 
the entanglement content of the ground state. 

PACS numbers: 03.67.-a, 03.67.Lx ,73.43.Nq 



I. INTRODUCTION 

Adiabatic Quantum Computation (AQC) has been 
proposed in 2001 by Farhi et al. ^ as an alternative 
way to perform quantum computation. It is based on 
the quantum adiabatic theorem which states that a 
system initially in its ground state, and subject to a suf- 
ficiently slow varying time dependent Hamiltonian will 
remain in the ground state as long as there is no energy 
level crossing in the course of the evolution. The protocol 
proposed by the authors consists in passing from a initial 
Hamiltonian Hq with a well known and easy to prepare 
ground state to the ground state of a final Hamiltonian 
(or problem Hamiltonian) Hp which encodes the answer 
to a given classical problem. This approach has been 
proved to be equivalent to the standard quantum gate 
model of computation The protocol of adiabatic 

computation amounts therefore to a controlled path in 
the space of hermitian operators acting on an Hilbert 
space of a physical system {H{s), s G [0,1]), which is 
the tensor product of n two level system (quantum bits). 
The path is usually taken to be the linear interpolation 
between Hq and Hp , resulting in a total Hamiltonian: 



H{s) = (1 - s)Ho + sHp. 



(1) 



The time T taken to perform the computation is such 
that the adiabatic theorem applies, with the probability 
of passing to an excited state remaining limited , which 
roughly translates to: 



T » A: 



(2) 



where A„ii„ is the minimum value of the energy differ- 
ence between the ground and the first excited states taken 
along the evolution. Therefore the time scaling with n 
(and so the computational efficiency) will be mainly de- 
termined by the behavior of the energy gap between the 
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two lowest energy states, and is rather long whenever 
this scaling is exponentially small. Up to now only ex- 
ponential decreasing gaps were proved to exist for some 
adiabatic protocols trying to solve NP-complete problems 
(BSHlEli) and some results, namely for the 3-SAT 
problem are rather inconclusive (P,|l[l3). For the un- 
sorted data base search a AQC algorithm was found that 
reproduce the gain of the Grover's algorithm for standard 
quantum computation pdl | . Even if no new improvements 
in algorithm design were obtained up to now using this 
approach, an advantage pointed out by is the robust- 
ness of this protocol against quantum errors. 



II. THE ADIABATIC ALGORITHM 

A. General Remarks 

Let us first discuss the general properties of the 
Hamiltonians considered in this paper. We scale all 
Hamiltonians such that they have a bounded spectrum. 
Even if for a typical physical system one expects the 
energy to be proportional to n, this is just a linear 
scaling and can be later taken into account in the total 
computational time for a realistic physical system. 

Given a classical computational problem one should 
first map its solution(s) onto a ground state of an 
Hamiltonian denoted Hp. We assume that the final 
measurement (output of the computation) is performed 
in the computational basis, and so, that the final Hamil- 
tonian is diagonal in this basis. Once Hp is defined, 
one can ask for an optimal initial Hamiltonian Hq and 
an optimal path H{s) in the parameters space (that 
maximizes the energy gap). 

In this paper we will only consider Hamiltonian 
paths of the simple form JQ. This choice is motivated 
by the following arguments. First it is clear that 
any Hamiltonian path of interest will be such that 
[dsH{s),H{s)] ^ 0. Indeed suppose that for a path 
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H{s) there exist an interval of values of s G [so,si] 
such that we have [dsH{s), H{s)] = 0, implying that 
Hamiltonians within that interval commute with each 
others. In this situation two cases may occur: either 
there is one or more level crossings between the ground 
state and the first excited state and, in that case, the 
adiabatic condition is no longer valid; or the ground 
state experiences no level crossing (so its correspondent 
eigenvectors remains the same) and the evolution 
amounts to a rescaling of the energies which could 
be performed "instantaneously" without breaking the 
adiabatic condition, because all the non-diagonal matrix 
elements of dsH{s) vanish. In this case the evolution 
between so and si needs not to be done adiabatically. 
In the following we suppose that all Hamiltonian paths 
do not have such "commuting" intervals and so we end 
up with a path such that [dsH{s), H{s)] ^ 0, V.s G [0, 1], 
which is the nontrivial part of the protocol. 

It should be clear that in most interesting cases 
the system will undergo a Quantum Phase Transitions 
(QPT) along the Hamiltonian path. Indeed the com- 
plexity of a classical computational problem translates 
into an increasing of the computational time with n. In 
the framework of the adiabatic theorem, this implies 
that there exists at least one value of s such that a 
gap closes (with increasing n) during the adiabatic 
evolution. Passing a vanishing gap corresponds to QPT 
in physical terms. In the vicinity of a QPT the system 
is described by its universality class which depends on 
the relevant couplings of the Hamiltonian. The family of 
Hamiltonians with the same relevant couplings presents 
the same behavior at the QPT and so share the same 
complexity when considered as adiabatic algorithms. 

The relation with QPT suggests that a given algo- 
rithm (path) should not fundamentally depend on 
" small details" of this path but rather on some relevant 
features near the QPT. This is why the simple interpola- 
tion Hamiltonian path can be chosen. Given Hp the 
choice of the path is then reduced to the choice of 

B. The Algorithm 

The output of the adiabatic computation is the ground 
state of the final Hamiltonian Hp. We build such a 
Hamiltonian by attributing to each possible classical con- 
figuration Xi, between the iV = 2" possible ones, a real 
value £i (energy) that measures how well the problem is 
satisfied by the string of bits xi. If solution of 

the problem we set e,; = 0, otherwise takes a non zero 
value, usually based on the problem statement (for ex- 
ample for the 3-SAT problem, is the number of clauses 
violated by the string Xi). The problem Hamiltonian 
reads 

H{s = \)=Hp=Y,e,\i)i^\, (3) 



where \i) runs over the iV = 2" states in the computa- 
tion basis which we take to be the tensor product of the 
individual eigenstates of cr^ for a two level system. If 
the physical system remains in its ground state, the final 
output corresponds to the minimization of the energy Si 
as a function of i. For sake of simplicity we take Hq to 
be diagonal in the x basis. Concerning the s parameter, 
we assume through this paper that it is a simple afiine 
function of the time such that s(0) = and s(T) = 1. 
Note that some authors have proposed to speed up the 
adiabatic evolution far from the QPT point leading to an 
efficient gain for the total computational time 

III. THE MODEL 

In the previously published papers on AQC two dif- 
ferent prescriptions of the initial and final Hamiltonians 
were used. The most common is an additive Hamiltonian 
made of interaction terms involving few qubits (pairs and 
triplets). Indeed this type of problems is usually given 
by a set of local constrains (concerning few variables e.g. 
3 for the 3-SAT problem). This additive prescription 
can also be used for iJp providing that in that case the 
ground state can be easily prepared. The other typ e 
of Hamiltonians are projector-like (H^ ~ i/) 0, IMEjI- 
A Hp of this type corresponds to an oracle-based 
problem which has two possible values of the energy: a 
(possibly degenerate) ground state energy whose states 
are the problem solutions and an excited energy for non 
solutions. 

We are interested in studying different types of adiabatic 
evolutions which differ in terms of the gap scaling. For 
that purpose we chose a fixed initial Hamiltonian Hq 
and we study a range of Hp with different gap scaling 
properties. Since H^ should be easy to implement and 
diagonal in the x direction, the most natural choice is a 
simple (normalized) magnetic field interaction along the 
X direction: 

^0 = ^ - -S.. (4) 
2 n 

where Sx = \ ■ We aim to present a simplified 

model of adiabatic computation by using an Hamiltonian 
H{s) which is solvable, while displaying some of the fea- 
tures which will eventually be found in the more realistic 
case, namely the spectrum diagonal in the z basis for 
_ff(s ~ 1), and a quantum phase transition at some in- 
termediate value Sc. We choose a hermitian operator 
with k spin interaction terms of the form of the following 
tensor product 

K = + (5) 

where i denotes the qubit. It is clearly a sum of fc-spin 
interaction terms (fc ranging from 1 to n) , whose strength 
depends on p. Introducing the total spin Sz = \ 'Yi!i=\ 
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and using the identity e""^^ = cosh(a)+sinh(a)crz, hp can 
be written as 



hp = (1 - p ) =6 



2\--| 2tanh"^(p)S5 



(6) 



Finally, introducing a ~ n tanli ^ (p) , we rescale hp into 
i/p in the following form 



gal _ g2f 



2 sinh (a) 



(7) 



Note that the n factor in the definition of a is intro- 
duced to obtain an Hamiltonian Hp such that nHp is 
an extensive operator. 

The ground state of this Hamiltonian is = 
which has zero energy. The mth excited states corre- 
spond to a state with m Is having a binomial degeneracy 
7 — "\| I for any finite value of a and whose energy 

[n—my.n: ^ 

depends on a (see below). For a —^ oo the Hamiltonian 
is proportional to a projection operator: 



Hp{oo) = (I - 



(8) 



and all the excited states have energy equal to one. The 
final Hamiltonian obtained this way is the Grover-like un- 

sorted data base searching considered in m and [m 

The limit a — > decouples the qubits and H(s) can be 
written as a sum of independent single qubit Hamiltoni- 



H{s,a^ 0) 



Jfe) 



(fe) 



(9) 



k=l 



which corresponds to a trivially separable problem that 
can be solved by parallelizing single qubit problems. 

Note that, although the above Hp ground state 
has a particular simple form, the properties described 
below would apply to any Hamiltonian obtained from 
H{s) under unitary transformations. Suppose a given 
problem has a solution corresponding to the ground state 
\w{J)) = ®k=i\'^k), where Jk € {0,1}, and the same 
energies and degeneracies as the Hp described above. In 
this case by performing the unitary transformation: 



U 



(10) 



to H{s), Hq remains invariant and the ground state of 
Hp is changed to \w{J)). Since U is unitary the gap scal- 
ing nature and the energy spectrum stay invariant. So, 
upon studying the particular form of Hp ((TJ, we have ac- 
cess to the gap behavior of problems which can have the 
whole set of possible 2" ground states i.e. solutions. The 
advantage of Hp{a) as a representant of this class is that 
it is symmetric under qubit permutation and so it can be 

written as a function of the total spin Sa = 5 X)fe '^^'^^ 
for a = x,y, z. Moreover since Hq and Hp{a) commute 



with 5*^ = X^a "^a ^''^d the ground state of Hp has max- 
imal value for the total spin J = the whole evolution 
will take place in this maximal spin sector spanned by the 
Dicke basis: n, where the two quantum num- 

bers stand for the value of the total spin {j = ^) and 
the spin projection along the z axes. All the other spin 
sectors can be disregarded because they are not coupled 
by the adiabatic evolution 



A. Density of states 

As explained above the complexity of the adiabatic 
evolution is related to what happens near the QPT, it is 
nevertheless interesting to describe the spectral proper- 
ties of H{s) for the extremal values of s (equal to zero 
and one). The density of states as a function of the en- 
ergy for s = and s = 1 can be given analytically for 
large n. For s = 0, considering the binomial degeneracy 
of each level, the density of states writes: 



(11) 



where Ms{ij->)du} is the total number of levels between uj 
and LJ -\- duj. For the case s = 1 one starts by remarking 
that the energies ujm in the maximal spin sector are given 
by: 



^m{s = 1) = 



„a(2f:-l) 



2 sinh a 



(12) 



where m is the level labeling m ~ 0, n. For large n one 
can invert this relation and obtain the energy density for 
the maximal spin sector: 



1 



a (1 — 2a; + cotha) 



(13) 



To get the density of states as a function of the energy 
for the total spectrum one has to consider the binomial 
degeneracy of each level: 



Ns=i{'jj)du = 



2n e-2"(- 



-2i^ainhQ) 1^2 
Jo -2 I 



TT a (1 — 2ijj + cotha) 



-duj. (14) 



The behavior of the density of states at s = 1 as a func- 
tion of n and a is shown in Fig. ^ Remark that for 
a = one recovers the Gaussian centered a.t uj = ^ re- 
sulting from the binomial degeneracy of the separable 
problem. For increasing a, the density of states is more 
peaked toward the value lu = I which characterizes the 
projector-like Hamiltonian. 



IV. MEAN FIELD APPROACH 

Since each qubit interacts in an equivalent way with all 
the other qubits, we expect that, in the thermodynamic 
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FIG. 1: Density of states as a function of the energy for s = 1 
for different values of ct. The solid lines correspond n = 100 
and dashed lines to n = 300, plotted for a € {0, 1, 2, 3}. 



limit, a mean field approach will give access to the exact 
ground state energy and signal quantum phase transi- 
tions whenever they occur. This is done by injecting the 
separable ansatz state: 



.Ucos(-)|0)+Sin(-K^|l) 



(15) 



and minimizing the energy in order to determine the free 
parameters and (p. Doing so we can see that a first 
order quantum phase transition occurs for some value of 
s = Sc{a){Ki 2.598) providing that a > FigElshows 




0.25^'"^ 



of (5*2;) is related to the abrupt change of the ground state 
component with respect to the fully polarized states in 

the X {\ ^) = ®r=i (75(10) + |1)))) and z directions 

(l^) = ®Li|0)). 

V. NUMERICAL ANALYSIS 

A. Energy Spectrum 

As predicted by the mean field approach the value of 
the ground state energy tends to a maximum around s = 
Sc{a) and approaches the mean field ground state energy 
for increasing n. At the critical point there is an energy 
level anti-crossing and the gap vanishes as n increases. 
Fig. 121 shows the energy spectrum (in the symmetric 
sector) as a function of s for different values of the a 
parameter and n. 
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FIG. 3: Energy levels as a function of s computed for different 
values of n and a. The series I presents the energy levels for 
n = 10 and a £ {0,3,10} respectively for The 
series // presents the energy levels for n = 20, the values of 
a are labeled as before. 



FIG. 2: Mean value of = — plotted in the a — s plane. 



the mean value of the observable , 



which presents 



a discontinuity along the first order quantum phase tran- 
sition line in the a ~ s plane. This line ends with a 
second order point (a,s) — (^^, 5—^ — 7 — ^. _i ) 

2+3\/6e7 sinh(^j 

where the values of observables are continuous non ana- 
lytic functions of a and s. For a > a discontinuity 



B. Analysis of the Two Less Energetic States 

For the cases where a phase transition exists {a > ac = 

we now analyze the behavior of the two less ener- 
getic states, in particular their projections along | =>) 
and I ff). Fig^ presents a zoom of the ground state 
anti-crossing (at s = Sc) with the first excited state an- 
teceded by some anti-crossings between more energetic 
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states. As n increases this cascade of anti-crossings gets 
closer to the QPT point. FigO shows the projections 
of the ground state and the first excited state along 
the fully polarized states | =J>) and | 11). We observe 
four different regions limited by the values of s where 
I iV'o/i) P or I (fr IV'o/i) P change abruptly. In the 
region > s > si there are several level anti-crossings 
between excited states but they do not affect significantly 
the first two states of the spectrum (FigQ)); we have 
\,Po) ^ I ^) (Fig© and IV'i) If,^). up to a very 
good approximation. At Si the first excited state suf- 
fers an anti-crossing with the second one and its projec- 
tion along I fl") increases drastically but remains different 
from one. FigElshows the projections onto the Dicke ba- 
sis of the ground state immediately after the anti-crossing 
and of the first excited state immediately before the anti- 
crossing (s ~ Sc). At s = Sc the transfer of components 
between the ground state and the first excited state is 
clearly manifested. For s > Sc the ground state increases 
slowly its projection along | =>) and at S2 the first excited 
state experiments another anti-crossing (Fig^, increas- 
ing drastically his projection along \j = = ^ — 1). 
For s > S2 there are no anti-crossings and the states in- 
crease slowly their projection along the ground and first- 
excited states: | ff) and \j = ^, m = ^ — I). 
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FIG. 4: Different regions where the ground state (blue) and 
the first excited state (red) undergo level anti-crossings, plot- 
ted for n = 30 and a = 5. 



VI. BEYOND THE MEAN FIELD 
APPROXIMATION 

In order to characterize further the system let us now 
look at the gap scaling (with n) near the QPT and the 
ground state entanglement content as measured by the 
concurrence [l^ . To study the quantum fluctuations 
around the mean field solution we can for instance ap- 
ply a method described in which uses the Holstein- 
Primakoff mapping from a given spin sector characterized 
by the value of S"^ and the algebra of boson operators, 
obtaining a interacting boson Hamiltonian which can be 
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FIG. 5: Square of the modulus of the inner product of the 
two lowest levels with the states | =>) and | ■fl) which are the 
ground state of H{s) for s = and s = 1 respectively. Here 
plotted for n = 50 and a = 5. (a) and (b) present respectively 
the evolution of the ground and first excited states projections 
along the Hamiltonian path, (c) combines (a) and (b), and 
clearly displays the exchange between these two states at s = 

Sc. 
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FIG. 6: Square modulus of the projections of the ground state 
immediately after the anti-crossing and of the first excited 
state immediately before the anti-crossing along the Dicke 
basis, plotted for n = 100 and a = 5. Note that the a loga- 
rithmic scale is used. 



expanded in powers of Using this method one is able 
to compute the ground state concurrence in the a — s 
plane and to predict the gap behavior based on the uni- 
versality class of the model. 
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A. Concurrence 

We measure the entanglement contend of the ground 
state by computing the concurrence of the (n — 2)-qubit 
traced density matrix. The entanglement is encoded in 
the finite size corrections so the quantity to study is the 
rescaled concurrence Cn = (TV — 1)C which is usually 
non trivial in the thermodynamic limit near a QPT. This 
quantity can be computed, for the symmetric spin sec- 
tor, as a function of the mean values of the spin operators 
i = x,y,z, In the present case the real na- 

ture of the density matrix leads to the simple expression: 

Cn = 1 : see [T^ ■ Fig[7| displays the concur- 
rence computed in the a — s plane. Note that along the 
first order line the concurrence is discontinuous in s. At 
the second order transition point this quant ity p resents a 
cusp like form as in the simple LGM model [13 ■ The en- 
tanglement entropy can also be computed for this model 
using the method discussed in and also displays a 
discontinuity along the first order line [T9l |. 




FIG. 7: Reduced concurrence in the thermodynamic limit ob- 
tained, by tracing over n — 2 arbitrary qubits, as a function of 
a and s. For « > «c the concurrence presents a discontinuity 
which increases with a. 



B. Scaling of the Gap 

Fig|51 shows the scaling of the minimun energy gap, 
obtained along the path H{s), with 71. For a < the 
gap is proportional to — . This arises because of our nor- 
malization choice of the total Hamiltonian. Had we cho- 



sen a normalization in which the energy was a extensive 
quantity (which can be obtained by multiplying the total 
Hamiltonian by n) and the gap would be constant in the 
large n limit. This result is simply derived from the stan- 
dard Holstein-Primakoff method. At the second order 
transition point {a — the numerical computation 

displays a clear slope crossover toward a still polynomial 
but non trivial exponent A ~ n~'^ (FigEJ- Numerically 
we find v to be close to | , which is the value found previ- 
ously in • for an also fully connected spin system but 
where the Hamiltonian was limited to two-body interac- 
tions. In the region where the first order QPT occurs 
a > the gap vanishes exponentially with n. This is 
a general behavior for first order QPT. 
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FIG. 8: Scaling of the energy gap at the QPT point with 
71 for different values of the a parameter. The curves range 
from a = to a = 3. One observes a clear crossover at the 
critical value Oc — (yellow) with a u value close to |. 



VII. CONCLUSION 

In this paper we discuss in detail an Hamiltonian 
evolution which should be viewed as a toy model for 
adiabatic computation. Indeed the phenomcnological 
properties of this system correspond to what is usually 
expected in more realistic implementations: a Hamilto- 
nian based on spin-spin interactions, a final Hamiltonian 
(s = 1) diagonal in the computational basis and a 
non trivial behavior for some intermediate values of s 
corresponding to a QPT. An interesting feature of our 
model is that it is build on a two dimensional parameter 
space which allows to trigger the Hamiltonian path 
form a trivial one (without QPT) to a regime where 
a first order QPT occurs, separated by a second order 
phase transition. The above phase space may serve as 
a template for more realistic cases. Given a classical 
computation problem the precise Hp formulation and 
the choice of the initial Hamiltonian Hq should result 
from an optimization process. A qualitative knowledge 
of the phase space is required for that analysis (in terms 
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of the topology of the first and second order phase 
transition manifolds). In particular an heuristic point 
of view would lead to looking (in the phase space) for 
second order QPT to built an optimum Hamiltonian 
path. 
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